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Abstract 

Convergence  of  linear  quadratic  regulator  (LQR)  problems  in  structures  is  discussed.  The  abstract 
formulation  of  the  system  using  a  variational  framework  based  on  sesquilinear  forms  is  considered.  Since 
convergence  theorems  require  uniform  stabilizability  of  the  finite  dimensional  approximating  system,  we 
present  a  detailed  proof  of  a  fundamental  lemma  due  to  Banks  and  Ito  which  can  be  used  to  easily  verify 
this  condition  for  many  applications.  Existing  results  for  the  well  posedness  of  the  infinite  dimensional 
system  and  convergence  of  Galerkin  approximations  are  summarized. 
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1  Introduction 

In  this  paper  we  discuss  in  detail  the  proof  of  Lemma  6.2  in  the  electronic  and  CRSC  tech¬ 
nical  report  versions  of  [1]  (stated  as  Lemma  7.13  in  [2])  which  allows  verification  of  uniform 
stabilizability  of  a  family  of  finite  dimensional  approximating  systems  arising  in  feedback  con¬ 
trol  formulations.  This  uniform  stabilizability  condition  is  sufficient  for  the  desired  convergence 
theorems  for  gains,  controls  and  trajectories.  Since  structural  applications  are  of  interest  to  us, 
the  partial  differential  equations  we  consider  are  second  order  in  time  and  are  motivated  by  a 
simple  example  of  a  cantilever  beam.  We  discuss  details  of  this  proof  since  it  has  appeared  only 
in  the  electronic  version  of  [1]  and  the  proof  appearing  there  was  only  a  sketch  of  the  arguments. 
We  believe  the  ideas  behind  this  proof  can  be  extended  to  treat  a  much  larger  class  of  examples 
than  those  indicated  below. 

The  system  we  consider  entails  Kelvin- Voigt  or  strong  damping  and  the  lemma  requires  a 
strong  assumption  regarding  the  relationship  between  the  stiffness,  Kelvin- Voigt  damping  and  air 
damping  coefficients.  This  assumption  is  reasonable  for  homogeneous  structures  such  as  beams, 
plates  or  shells  with  actuators  (such  as  piezoceramic  patches)  embedded  in  a  manner  so  that 
models  with  material  properties  which  do  not  vary  across  the  region  of  the  actuators  are  good 
approximations. 

The  theoretical  control  results  we  present  are  valid  for  systems  with  bounded  observation 
operators.  Since  current  sensing  devices  yield  observation  operators  that  are  discontinuous  in 
nature,  we  refer  the  reader  to  extended  results  regarding  unbounded  observation  operators  in  [1, 

1Research  supported  in  part  by  the  U.S.  Air  Force  Office  of  Scientific  Research  under  grants  AFSOR  F49620- 
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Section  6.2].  We  also  assume  no  persistent  exogenous  force  is  acting  on  the  system  and  thus  we 
are  interested  in  applications  in  which  the  structure  or  actuators  start  with  a  deformation  and 
vibrates  to  a  steady  state.  This  is  also  useful  in  applications  where  the  system  is  subjected  to 
an  impulsive  force.  For  systems  with  exogenous  disturbance,  the  theory  is  less  complete  but  a 
number  of  computational  investigations  have  been  carried  out  on  systems  with  periodic  exogenous 
forces.  We  refer  the  reader  to  [2]  for  a  summary  of  results  for  infinite  and  finite  dimensional 
control  of  systems  with  periodic  exogenous  disturbance.  Numerical  confirmation  involving  LQR 
control  of  thin  cylindrical  shells  with  transient  and  periodic  exogenous  disturbances  can  be  found 
in  [3,  4,  5].  For  experimental  and  numerical  results  of  control  of  plate  systems,  see  [6]. 

In  Section  2  we  discuss  the  abstract  system  and  introduce  the  motivating  example  which  is 
control  of  transverse  vibrations  of  a  beam  as  presented  in  [2].  We  then  summarize  results  for 
infinite  dimensional  control  in  Section  3.  Approximation,  LQR  control  and  the  detailed  proof  of 
Lemma  6.2  of  [1]  are  given  in  Section  4. 


2  Abstract  System 


Consider  the  abstract  second-order  (in  time)  variational  system 

(w(t),'ip)v,v  +  a2(w(t),‘ip)  +  ay  (w(t),ip)  =  ( Bu{t),ip)v.'V 
w(0)  =  wq,  w(0)  =  w\  , 


where  oq  and  a2  are  sesquilinear  forms  from  R  x  R  to  C.  Let  R  and  H  be  complex  Hilbert 
spaces  with  R  continuously  embedded  in  H  (i.e.,  U\\h  <  k\\<t>\\v),  forming  a  Gelfand  triple 
R  <->  H  <->  R*  (for  details  regarding  Gelfand  triples,  see  [7,  8]).  Here  R*  and  H*  are  the  dual 
spaces  to  R  and  H ,  respectively  and  H  is  identified  with  H*  through  the  Riesz  map.  We  take 
the  duality  product  (•,  -)y*xv  on  R*  x  R  to  be  the  unique  extension  by  continuity  of  the  scalar 
product  of  H  defined  on  H  x  R.  Thus,  the  elements  v*  £  V*  have  the  representation 

v*(v )  =  {v*,v)v,  v.  Furthermore,  assume  that  the  embedding  i  from  R  into  H  is  compact,  and 
that  the  stiffness  sesquilinear  form  oq  is  R-continuous,  positive  and  symmetric,  i.e., 


(HI) 

(H2) 

(H3) 


Wii<t>,ip)\  <  Cl  \<t>\y  iV’ly,  for  some  a  €  IR 
Reoi{4>,(j>)  >  c2  \cj>\y,  for  some  c2  >  0 

o-i  (< t>A)  =  bP,4>) 


(R-continuous) 
(R-elliptic) 
(symmetric)  . 


Assume  further  that  the  damping  form  a2  has  the  properties 

(H4)  \a2{4>,'ip)\  <  cs  \<j)\v  IV’ly,  for  some  C3  €  IR  (R-continuous) 

(H5)  Rea2((f),(j))  >  C4  IV’ly,  for  some  C4  >  0  (R-elliptic)  . 

From  the  continuity  properties  (HI)  and  (H4),  we  obtain  operators  Ai,A2  £  £(R,  R*)  defined 
by 


(Al<j),  1p)y.  y  =  (7l{(t>,.i>)  >  Vl),^R 

(' A2<f>,1p)y,ty  =  a2(<t>,ll>)  ,  VV>  ,V»GR, 

and  hence  an  equivalent  formulation  of  (2.1)  is  given  by 


(2.2) 


w(t)  +  A2w(t )  +  A\w(t)  =  Bu(t)  in  V* 
w(0)  =  wq,  w(Q)  =  w\  . 


(2.3) 
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Figure  1:  Cantilever  beam  with  piezoceramic  patch  pair. 


The  control  operator  B  :  U  — »  V*  where  U  is  the  control  input  space  is  typically  unbounded  due 
to  the  discontinuous  geometry  of  actuators.  The  input  u(t)  usually  models  voltage  input  to 
the  actuators  in  smart  material  structure  applications. 

The  system  (2.1)  or  equivalently  (2.3)  arises  in  the  abstract  formulation  of  partial  differential 
equations  governing  smart  material  structures.  To  illustrate,  consider  the  transverse  vibrations 
of  a  homogeneous  beam  with  length  £,  thickness  h,  width  6,  linear  mass  density  p,  Young’s 
modulus  E  and  Kelvin  Voigt  damping  cd ■  Assume  that  a  pair  of  identical  piezoceramic  patches 
are  bonded  to  opposite  sides  of  the  beam  covering  the  region  y\  <  y  <  y2  (see  Figure  1).  We 
denote  the  Young’s  modulus,  linear  mass  density,  thickness  and  damping  coefficient  of  the  patches 
by  Epe, ppe,hpe  and  cd  ,  respectively.  Cantilever  end  conditions  are  assumed  with  the  fixed  end 
at  y  =  0  and  free  end  at  y  =  l.  Transverse  displacements  of  the  beam  will  be  denoted  by  w 
and  air  damping  coefficient  by  ca.  Force  and  moment  balancing  yields  the  strong  form  of  the 
equation  (see  [2]  for  details  in  the  derivation)  which  when  written  in  weak  form  is  given  by 


fl  \~d2w  ~d2wd2<f>  - — -  d3w  d24> 

J0  \PW<t>  +  EI  dy2  dy 2  +  CnIdy2dt  dy2 


dw  ,  . 

+  Ca-g^4>-  (bMv)pe 


&±\ 

dy 2  J 


dy-  0  , 


for  all  4>  €  H'K 0,£),  where 


(2.4) 


Hl( 0,i)  =  H2(0JM0)  =  0'(O)  =  0}  . 


(2.5) 


Due  to  the  presence  of  the  patches,  the  linear  mass  density  p(y)  =  phb  +  2bppehpexPe(y)  is 
piecewise  constant  with  the  characteristic  function  xPe(y)  used  to  isolate  patch  contributions. 
Here 

EI(y)  = 
cDl{y)  = 

and  a3  =  (h/2  +  hpe)3  —  h3 / 8. 

The  external  moment  ( bMy)pe  depends  on  the  voltages  supplied  to  the  two  patches.  Denoting 
the  outer  and  inner  patch  voltages  by  V\  ( t )  and  V3  ( t ) ,  respectively,  the  external  moment  is  given 
by 

(bMy)pe  =  -KBXpe{y)  [Vi (f)  -  V2(t)\  ,  (2.7) 


K3b  2 b„  .  , 

^  12  3  ^pea3Xpe\y) 

h3b  2b 

°d~[2  +  ~JCDpea3Xpe^y' 


(2.6) 
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where  KB  =  EpebdS\(h  +  hpe)  depends  on  the  piezoceramic  material  properties.  Coupled  to 
the  system  are  the  cantilever  boundary  conditions 

w(t,0)  =  ^(t,0)  =  0,  My(t,t)  =  ^ My(t,t )  =  0  (2.8) 

and  initial  conditions 

w(0,y)  =  w0(y),  -^(0, y)  =  wx{y)  .  (2.9) 

Here  My  is  the  internal  moment  resultant  with  expression 


My 


1 

b 


~d2,w  - — -  d2w 

EI~dy2  +  °DlWdi 


To  abstractly  formulate  the  weak  form  (2.4),  let  V  =  i?|(0,  t)  and  H  =  L2(0,l).  It  readily 
follows  from  standard  Sobolev  theory  that  V  is  continuously,  densely  and  compactly  embedded 
in  H  and  hence  with  H  forms  a  Gelfand  triple  V  c— ^  H  c->  V*.  For  ij>,  ip  £  H,  define  the  H  inner 
product  to  be 


(<P 


*I>)h  =  / 

Jo 


phb(f>ipdy  , 


and  for  <j),ip  £  V,  the  stiffness  and  damping  sesquilinear  forms  are  defined  by 


o-i  {<t>,ip) 
(<t>,1p) 


d2cf>  92tl> 
dy 2  dy 2  dV 


f 

J  0 


cDi(y) 


d2cf>  d2ip 

Thflhf 


dy  +  ca 


(2.10) 


The  control  operator  B  :  U  V*  is  given  by 

fe  d2ib 

(Bu(t),t[>)v,  v  =  (bMy)pe -Qy^dy  .  (2.11) 

Here  U  =  IR  and  the  vector  u(t)  £U  represents  the  time  varying  voltage  to  the  inner  and  outer 
patches.  Finally,  we  define  the  l?£-equivalent  inner  product  on  V  using  the  stiffness  sesquilinear 
form,  i.e.,  (■,  -)v  =  <7i (■,■)■ 

It  can  be  easily  shown  that  the  sesquilinear  forms  (2.10)  satisfy  (H1)-(H5)  and  thus  the 
weak  form  (2.4)  with  corresponding  boundary  (2.8)  and  initial  conditions  (2.9)  can  be  abstractly 
formulated  using  (2.1)  (equivalently  (2.14)). 

To  obtain  a  first  order  formulation  of  the  system  (2.1)  (or  equivalently  (2.3))  amenable  to 
semigroup  analysis  and  control  methodologies,  we  define  the  product  spaces  T~L  =  V  x  H  and 
V  =  V  x  V  with  norms 

||(0i,02)||^  =  1101  lly  T  ||02 ll^f 

11(01, 02 ) 1 1 v  =  1101  lly  +  1102 II v"  ■ 

It  can  be  readily  verified  that  these  product  spaces  also  form  a  Gelfand  triple  V  4  H  4  V”, 
where  V*  =  V  x  V* .  The  control  operator  is  then  reformulated  as 


Bu(t)  = 


0 

Bu(t ) 


and  we  define  a  sesquilinear  form  a  :  V  x  V  -4  C  by 


<j($,  ^)  =  -  (02 , 01 )  V  +  <71  (01, 02 )  +02(02,02) 
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for  $  =  (4>i,4>2),  4/  =  (^1,^2)  €  V.  Thus,  for  z(t )  =  (w(t),w(t)),  the  second-order  system  (2.1) 
formulated  in  first-order  form  is 


(z(t),  *>v.,v  +  *)  =  ^>v*,v  .  VfeV 

2(0)  =  z0  =  (w0,wi)  . 


(2.12) 


The  V-continuity  of  a  and  V-ellipticity  of  <r(-,  ■)  +  A{-,  -)n  for  A  >  0  was  detailed  in  [2,  p.109]. 
This  guarantees  the  existence  of  the  operator  A  £  C(V,  V*)  given  by 


a($,$)  =  <^,$>v.iV  • 


To  write  (2.12)  in  an  equivalent  strong  form,  we  restrict  A  to  the  system  operator 


domA  —  {(<^1,^2)  £  £  V,Ai<pi  +  A^fa  G  i?} 


(2.13) 


where  Ai,A2  are  defined  in  (2.2).  It  should  be  noted  that  A  is  the  negative  of  the  restriction  to 
domA  of  the  operator  A  so  that  ct($,  \H)  =  (-Al$,  4/ )n  for  $  £  dom^4,  ft  e  V.  A  strong  form  of 
the  abstract  system  model  (2.12)  is  given  by 


z(t)  =  Az(t)  +  Bu(t)  in  V* 
z( 0)  =  z0  • 


(2.14) 


Existence,  uniqueness  and  continuous  dependence  on  data  of  the  solution  to  (2.1)  was  first 
proven  in  [9]  and  can  also  be  found  in  [2,  Chapter  4].  For  the  first  order  form  of  the  system, 
existence  and  uniqueness  of  the  solution  together  with  its  equivalence  to  the  second  order  weak 
solution  are  presented  in  [2,  9]. 


3  Infinite  Dimensional  Control 

The  abstract  model  (2.12)  (equivalently  (2.14))  is  useful  in  applications  where  the  structure 
undergoes  initial  deformations  and  vibrates  to  a  steady  state.  Control  methods  for  this  case  are 
designed  to  attenuate  only  the  transient  state  responses. 

The  output  observations  in  the  observation  space  y  are  given  by  z0b  =  Cz(t)  where  C  £ 
lz{T~L,y)  is  bounded.  As  already  noted,  the  operator  C  is  often  unbounded  in  applications  but 
we  only  consider  the  bounded  case  here.  We  also  make  the  simplifying  assumption  that  the 
full  state  z  =  (■ w,w )  is  available  for  the  computation  of  the  feedback  control  u(t).  In  many 
practical  applications,  current  measuring  devices  can  only  deliver  partial  state  measurements 
hence  compensators  must  be  included  in  the  control  design.  The  reader  is  referred  to  [2,  Chapters 
7.5  and  8]  for  discussions  regarding  compensators  and  [1]  for  a  summary  of  results  on  unbounded 
observation  operators  C. 

The  quadratic  functional  we  minimize  in  order  to  determine  the  optimal  control  u  for  the 
infinite  horizon  control  problem  is 

J(u,zo)  =  J  j||C,z(f)||y  +  Hl/2u(t)  ^ |  dt  (3.15) 
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subject  to 


Az(t)  +  Bu(t ) 


z{t)  = 
z( 0)  =  z0  ■ 

Here,  the  positive,  self-adjoint  operator  1Z  =  (771/2)2  €  C{U,U)  is  used  to  soft  constrain  the 
control  input.  We  do  not  state  results  for  the  finite  horizon  problem  but  the  reader  is  referred 
to  discussions  in  [1,  Theorem  3.1]  and  [2,  Chapter  7.2.1]. 

We  first  give  the  definitions  for  the  pair  (A,B)  to  be  stabilizable  and  (A,C)  to  be  detectable 
before  stating  a  theorem  which  which  uses  these  conditions  to  guarantee  the  existence  of  optimal 
controls  minimizing  (3.15). 

Definition  3.1  The  pair  (^4,  B)  is  said  to  be  stabilizable  if  there  exists  an  operator  K  £  £(V* ,  U) 
such  that  A  —  BK.  generates  an  exponentially  stable  semigroup  on  V*  (i.e.,  there  exists  M  >  1,  u  > 
0  such  that  || e*!”4  e/c)  ||£(V„)  <  Me~wt). 

Definition  3.2  The  pair  ( A ,  C)  is  said  to  be  detectable  if  there  exists  an  operator  T  €  £((F,  V*) 
such  that  A  —  TC  generates  an  exponentially  stable  semigroup  on  V*. 

Theorem  3.1  If(A,B)  is  stabilizable  and  {A,C)  is  detectable,  then  the  algebraic  Riccati  equation 

{A*H  +  JlA-B.BTZ-1B*Jl  +  C*C)z  =  Q  VzeV  (3.16) 

has  a  unique  non-negative  solution  n  €  £(V*,V),  A  —  BTZ~lB*Tl  generates  an  exponentially 
stable  closed  loop  semigroup  S(t)  on  77,  V,V*,  and  the  optimal  control  that  minimizes  (3.15)  is 
given  by 

u(t)  =  -7 z-xB*m(t) 

where  z(t)  =  S(t)zo  for  zo  €  V*. 

Note:  Theorem  3.1  is  Theorem  7.5  in  [2]  and  its  proof  can  be  found  in  [1]. 

4  Approximation  and  Finite  Dimensional  Control 

The  solutions  to  (2.1)  or  (2.14)  together  with  the  optimal  controls  given  by  Theorem  3.1  are 
infinite  dimensional.  For  numerical  applications,  we  use  Galerkin  approximations  to  obtain 
solutions  in  finite  dimensional  subspaces  V ^  C  V  C  71.  The  bases  for  these  subspaces  can 
consist  of  modes,  splines,  polynomials,  finite  elements,  or  reduced  basis  elements. 

Since  V  =  V  x  V,  we  use  the  superscript  M  to  denote  W-dimensional  subspaces  of  V  and  use 
the  superscript  2W"  to  denote  the  2A/"-dimensional  subspace  V2A  =  V'Kr  x  V'Kr .  The  following 
approximation  condition  is  necessary  for  the  convergence  of  the  finite  dimensional  solutions  to 
the  infinite  dimensional  solutions  to  (2.1)  or  (2.14) 

(H1N)  Let  V  be  a  Hilbert  space.  For  any  zeV,  there  exists  a  sequence  z ^  €  \N  such  that 
\z  —  z*r jv  ->0asiV^oo. 

We  now  define  the  operators  A±r  :  V ^  >  V N  and  A %  :  V ^  -¥  V ^  which  approximate  A\ 
and  A‘2 ,  respectively,  by  restricting  the  corresponding  sesquilinear  forms  to  V ^  x  VAr ,  i.e., 

{Ai  <P,ip)H  =  ui  (<t>,i))  'i4>,%!)£VN 

(4.17) 

(A^4>,il))H  =  a2  V<j>,il>  &VN  . 
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Similarly,  the  operator  A2J v  :  V2A/”  — »  V2Ar  is  defined  by  the  restriction  of  a  on  V2"v  x  V2A/”  with 
definition 

(-A2M§,  V)n  =  ct($,  $)  V$,feV2Af.  (4.18) 

It  readily  follows  that 


(4.19) 


(4.20) 


2Af  __  0  I 

-A?  -A? 

The  usual  projection  operators  from  H  onto  V ^  and  from  %  onto  V2A/”  are  denoted  by  V ^  and 
V2jKf,  respectively,  and  are  defined  by 

VMf>  G  VM  and  (VN  -  <j>,  ip)H  =0  €  VM 

V2M§  g  V2AT  and  =  0  V4/  G  V2Ar  . 

The  control  operator  B  is  approximated  by  B2Ar  by  restricting  it  to  the  finite  dimensional  subspace 
using  its  adjoint 

() B2Nu ,  $)w  =  ,  V*  G  V2M  , 

and  C2jV  is  the  restriction  of  C  to  V2Ar .  Thus,  the  finite  dimensional  analogue  to  (2.14)  is  given 
by 

z2X(t)  =  A2Mz(t)+B2Mu{t)  mV2M 
z2M{ 0)  =  V2Xz0  • 

The  following  is  Lemma  4.1  in  [1]  which  guarantees  the  convergence  of  the  finite  dimensional 
Galerkin  solutions. 


(4.21) 


v*,v  ’ 


V*  G  V 


2A r 


(4.22) 


Theorem  4.1  (Lemma  4-1  in  [I])  Suppose  (H1N)  is  satisfied  and  let  Bu  G  L2(0,  T;  V*)  and 
zo  G  %.  If  z(t)  G  V  is  the  solution  to  (2.14)  and  z2J^(t)  G  V2jKr ,t  >  0  satisfies 

jt  (z2*(t),  *)n  +  a(z2*(t),  *)  =  (B2*u(t),  i>) 

z2N( 0)  =  V2Nz0  , 

then  the  error  function  e2jKf(t)  =  z2J^ (t)  —  z(t )  satisfies 

\\e2M(t)\\n^D 

and 


f  ||e2jV"(s)||y  ds  ->  0  as  N oo 

Jo 


uniformly  in  t  G  [0,T]. 


The  approximate  infinite  horizon  control  problem  for  the  system  involves  finding  the  control 
G  T2(0,  oo;  U)  which  minimizes 


2A ft,,  J2H 


J22V( 


subject  to  z2J v  satisfying 


z°M) = + ||k,/2»(*)|[} 

z2N(t)  =  A2MZ2M  +  B2Xu(t)  ,  t  >  0 

22AT(0)  =  V2NZ0=Z2N 


■  dt 


(4.23) 


(4.24) 


The  analogous  definition  of  stabilizable  and  detectable  matrix  pairs  for  finite  dimensional  systems 
will  be  stated  before  we  present  the  convergence  theorem. 
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Definition  4.1  The  pair  (A2J^  ,B2J^)  is  said  to  be  uniformly  stabilizable  if  there  exist  constants 
Mi  >  l.uJi  >  0  independent  of  N  and  a  sequence  of  operators  K?J v  £  jC(V2J^,U)  such  that 
supjy-  J/C2Ar|  <  oo  and 

et(A™ -B™  K?u)<plM  z 

for  Z  £Li. 

Definition  4.2  The  pair  (A2jKf  ,C2jKf)  is  said  to  be  uniformly  detectable  if  there  exist  constants 
M2  >  1, 0^2  >0  independent  of  N  and  a  sequence  of  operators  T21^  £  jC(y,V2Af)  such  that 
supjy-  \Jr2J^\  <  oo  and 

gt  (A2U -F2MC2M)<p2M  z 

for  z  £TL. 

Theorem  4.2  Suppose  V  H,  where  the  embedding  i  is  compact.  Let  the  sesquilinear  form 
a  associated  with  the  first-order  system  (2.14)  be  continuous  and  V -elliptic.  Assume  that  the 
operators  A,B,C  are  such  that  ( A,B )  is  stabilizable  and  ( A,C )  is  detectable  where  B  £  £(U,V*) 
is  unbounded  and  C  £  £{TL,  0^)  is  bounded.  Consider  an  approximation  method  which  satisfies 
( H1N ).  Finally,  suppose  that  for  fixed  No  and  N  >  No ,  the  pair  (A2jKf ,  B2Af)  is  uniformly 
stabilizable  and  (A2J^ ,C2jKf)  is  uniformly  detectable. 

Then  for  N  sufficiently  large,  there  exists  a  unique  nonnegative  self-adjoint  solution  II2Ar  £ 
jC(V*,V)  to  the  Nth  approximate  algebraic  Riccati  equation  (3.16)  in  V2A/”  with  A,B,C  replaced 
by  A2Ar ,B2^ ,C2J^ ,  respectively.  There  also  exist  constants  M3  >  1  and  u> 3  >  0  independent  of 
N  such,  that  S2Ar(t)  =  e(^  sat.sfies 

\\S2X(t)\\v2M  <  M3e-U3t  ,t>  0 

or  equivalently 

t  (  A2U -B2Un~ 1 B2**  ) 
e  v  zo 

Additionally,  the  convergences 

U2ATV2ATz  4  n z  in  V  for  every  z  £  V* 

\\B2*r*n2ArV2Ar  -  B*U\\c(n  u)  -s-  0  , 

as  N  — >  00,  of  the  Riccati  and  control  operators  are  obtained.  Moreover,  the  feedback  system 
operator  A  —  BIN1 B2^  II2Ar  generates  an  exponentially  stable  analytic  semigroup  onTL  and  for 
every  z0  £  %, 

J2N  (-B2M’U2N  z(-),z<^  -  J(u,z0)  <  e(N)  \\z0\\2n 
where  e(N)  -»  as  N  — >  00. 


<  M3e  W3t\\zo\\H  ,t>0,zo£'H. 


<  M2e  W2t  \\z\ 


m 


<  M\e 


-Wit  I 


m 


Theorem  4.2  is  Theorem  7.10  in  [2]  in  which  the  assumptions  of  uniform  stabilizability  of 
(A2^,  B21 v")  and  uniform  detectability  of  ( A 2A/”,  C2Af)  can  be  difficult  to  directly  verify.  Hence  we 
need  additional  results  based  on  readily  confirmed  assumptions.  For  general  first-order  systems, 
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Lemma  4.7  in  [1]  (which  is  restated  as  Lemma  7.12  in  [2])  guarantees  that  (A2Af ,  B2^)  is  stabi- 
lizable  provided  (A,B)  is  stabilizable  and  the  injection  V  %  is  compact.  For  the  second-order 
system  rewritten  as  a  first-order  system,  the  definition  of  the  product  spaces  V  =  V  x  V  and 
'H  =  H  x  V  precludes  the  possibility  of  compactness  of  i  :  V  <— >  R,  even  when  V  H  is 
compact.  Hence  we  require  the  following  lemma  to  obtain  uniform  exponential  bounds  on  the 
approximating  semigroups  S2Af(t). 

Lemma  4.1  (Lemma  6.2  in  [1])  Suppose  V  A  H,  where  i  is  compact.  Moreover,  suppose  that 
the  damping  sesquilinear  form  a2  can  he  decomposed  as  a2  =  7 ay  +  a2,  for  some  7  >  0,  where 
the  continuous  sesquilinear  form  a2  satisfies  for  some  p  £  IR 

Rea2(<t>,  4>)  >  --yMWv  ~  lMWh  for  all  (f>  £  V  . 

Finally,  suppose  that  the  operator  Tlf1  A2 ,  where  A2  £  C(V,V*)  is  defined  by 

(-4 <f>,  v)y,tV  =  ai{4>,  v)  (4.25) 

is  compact  on  V. 

Let  T  denote  the  open  loop  semigroup  generated  by  the  product  space  operator  A  and  let  T 2Af 
be  generated  by  A2Af .  If  for  some  v  >  p  and  M  >  1 

mt)\\c(U)  <  Mevt  ,  t  >  0  ,  (4.26) 

then  for  any  e  >  0  there  exists  an  integer  Me  such  that  for  M  >Me 

\\T^N(t)V2M\\cm  <  Me(i/+e)t  ,  t  >  0  , 

for  some  constant  M  >  0  independent  of  N. 

Since  the  proof  of  Lemma  6.2  in  [1]  is  given  only  in  sketchy  form,  we  give  a  detailed  proof  here. 
We  first  state  and  prove  the  following  lemmas  which  will  be  used  in  proving  Lemma  4.1  below. 

Lemma  4.2  (A^)  1  =  P^A^1  on  VAr ,  where  P :  V  ->  V ^  is  defined  by 

<7i  {P^<f>  -  <j),  ipN)  =0  for  all  £Vn  ,<t>£V  .  (4.27) 


Proof: 


•  Since  <7i(-,  ■)  gives  rise  to  an  inner  product  on  V  whose  induced  norm  is  equivalent  to 
the  usual  norm  in  V,  then  a  closer  inspection  of  (4.27)  reveals  that  P(f  is  the  projection 
operator  from  V  into  V ^  under  the  cry  (■,  ■)  inner  product,  i.e.,  the  projection  of  Vi  onto 
VAf  where  V±  is  V  with  the  oy  inner  product.  It  follows  immediately  that  P 'y  from  V  onto 
VAf  is  well  defined  and  linear.  Now  let  (j/^  be  any  arbitrary  element  of  VAr .  Then  for  any 
£  VN , 


(A^4>NAN) 


{<t>N^N)v.,v  *  (from  (2.2)) 

(< bM,ipAr)H  ,  (since  <f>^  £  V) 

(A? 

(Ji[{Ai)  1  <pAr,ipAr')  ,  (from  (4.17))  . 
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We  note  that  Ai  is  invertible  in  V ^  follows  from  (4.17)  and  the  P-ellipticity  of  o\ .  Thus, 

«ji  (Ailcj)M  -  (. A f )  1  < A,  A)  =  0 

for  every  ^/A  G  V'Kr .  But  by  the  defining  the  relationships  (4.27)  we  have 

<7i  (i^AA  -  A  AAA)  =  0  for  all  A,  A  G  VM  . 

It  follows  that  Pyf  A*  A  =  (AO  1  A  and  this  completes  the  proof.  ■ 

Lemma  4.3  (A")  1  A"  =  Pyf  A^1  A2  on  VM,  where  A f  is  defined  by  the  restriction  of  02  on 
VN,  i.e.,  Af  :  VAr  V  VN  is  defined  by 

(AA  4>)h  =  A  A  m  V<t>,  ^  G  PA  (4.28) 

Proof: 

•  Let  A  £  PA  Then  for  all  A  G  VN 

a  1  ^(AO  1  if  A ,  A)  =  ^if  A  >A)  y  (from  Lemma  4.2, 

(4.27)  and  (2.2)) 

=  (if  A,  A)  (since  if  A  G  VN) 

=  ^(^j  A)  (from  (4.28)) 

=  (from  (4.25)) 

\  ivy 

=  (aa1^^,/) 

\  Ivy 

=  <7i  A)  (from  (2.2))  . 

Therefore,  cy  ((AO  1  _  ,  A)  =  0  ,VA"  G  PA  We  then  conclude 

from  (4.27)  that  (AO _1  if  A  =  Pf  A1-^  A  VA  G  PA  ■ 

Lemma  4.4  ||  A*  -  Pf  A*  Il,c(y,y)  0 
Proof: 

•  First  we  consider  A 1  as  a  compact  operator  from  P  to  P  by  using  the  compact  injections 
i  :  V  <->  PI  and  i*  :  H*  <->  P*  and  setting  A*  =  A^Pi.  Next,  note  that  the  convergence 
Py[w  — >  Pc  ,  Vu>  G  P  is  evident  from  (H1N)  where  Pf  is,  as  introduced  above,  the  pro¬ 
jection  operator  from  V\  onto  VAr  (again  Pi  is  P  with  the  the  norm  induced  by  the  <Ji(-,  •) 
inner  product).  Thus  convergence  follows  from  Pf  z  —  z\  <  zAf  —  z |  ^  <  fc  J  A  —  z\y. 
Then 

V11  ~  PvA^Wcyv)  =  sup  (Ai1  -  P^AiAv 

1 '  1  Nk<i  V  J  v 

=  sup  1(7-  PfMAIIy 
\\v\\v<l 

=  sup  \\{I  -  P^)w\\ 

w€U 
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where  U  =  A1 1  ^{||u||y  <  1}^  is  a  relatively  compact  subset  in  V  since  A1 1  is  a  compact 
operator  from  V  to  V.  Thus,  by  Chatelin’s  Theorem4, 

sup  \\(I  -  P%)w\\v  ->0  . 

WSzU 

This  gives  us  the  desired  result 

II A  1  -  KA  1||/:(viv)  0  ' 

■ 

Lemma  4.5  A,  1 A  -  P{fA,  1 A2  — >  0  where  Ai  1 A2  is  compact  on  V . 

Vl  c(v,v) 

Proof: 

•  The  proof  of  this  lemma  is  similar  to  the  proof  given  for  the  previous  lemma.  ■ 

We  now  give  the  proof  of  Lemma  4.1. 

Proof  of  Lemma  4.1: 

•  We  first  express  the  sesquilinear  form  a  in  terms  of  components, 

<j($,  $)  =  -  (fa,'lpl)Vl  +Vl(<t>l,i>2)  +V2(<h,1p2)  , 

then 

Rea($,$)  =  Re{-{4>2,4>i)Vl+<Ti(4>u4>2) +02(4>2,<f>2)} 

=  Reu2 (<^2 ,  $2)  ,  since  is  equivalent  to  a i(-,-) 

=  yiJeui  ((j)2,  $2)  +  Rea2(cp2,  <h)  >  since 

<72(1/7  ip)  =  7<7i  (ip,  ip)  +  <72  (Vb  ip)  for  some  7  >  0 

>  ic2\\cp2\\v-  ^iWhWv  -  tA\<PA2H  >  since 

Reai(<p,<p)  >  c2||0||y  (see  (H2)),  and 
Rea2(ip,ip)  >  -^j-Wipllv  ~  lA\i>\\2H 

=  ~  V\\fo\\H 

=  ‘jiWfaWv  +  »Mi\\h  -»{\\<h\\2H  +  \\<Pi\\h} 

>  min{^^}  ll$llv -»\\®\\h  ■ 

It  follows  that  the  linear  operator  A  associated  with  <7  (see  (2.13))  is  the  infinitesimal 
generator  of  a  Co  semigroup  T (t)  satisfying 

WmWcm  <  Me"*  ,  (4.29) 

4Theorem  3.2  in  [10]:  Suppose  that  Tnx  4Ti,i£X.  Then  for  any  relatively  compact  set  U,  sup^gy  I \{Tn  ~ 
T)x ||  ->■  0. 
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and  thus  (4.26)  holds  with  v  =  p.  Now  suppose  (4.26)  holds  for  any  v  >  /i,  i.e.,  suppose 
||T(f)|| an)  —  Mevt  is  satisfied.  To  prove  the  lemma,  we  show  Me  >  0,  there  exists  Ne  such 
that  for  N  >  Ne 

||TaA/’(t)PaW'||  <  Me{v+£)t  ,  t  >  0  .  (4.30) 

The  next  step  involves  writing  the  semigroup  T(t)  in  terms  of  the  generator  A  using  the 
inverse  of  the  Laplace  transform.  From  the  result  above  that  A  is  the  infinitesimal  generator 
of  a  Co  semigroup  satisfying  (4.26),  it  follows  that  A  —  vl  generates  a  Co  semigroup  of 
contractions  S(t)  satisfying  ||<S(t)||;C^  <  M,  where  T{t)  =  S(t)evt .  Theorem  6. A  in  [7] 
guarantees  the  existence  of  6  such  that  0  <  S  <  n/4  and 

p(A  -  vl)  D  Ej  =  {z  €  (C  :  |arg(z)|  <  7r/2  +  <5}  . 

By  Theorem  7.7  in  [11],  we  can  express  S(t)  as 

S(t)  =  -5-t  [  ext(\I-  {A-vI))-ld\ 

2 7TI  Jy 

where  T  is  a  smooth  curve  in  E5  U  {0}  running  from  ooe~*’3  to  ocel§  for  7t/2  <  ()  <  -k/2  +  5. 
Since  T(t)  =  evtS(t ),  we  have 

T(t)  =  — — t  [  evtext{{\  +  v)I  -  A)-ld\  . 

2m  Jy 

Shifting  the  path  of  integration  T  by  v  and  denoting  it  by  V,  we  have 

T(t)  =  J~[  ex\\ I  -  A)-ld\  .  (4.31) 

2m  Jy> 

The  finite  dimensional  semigroup  can  similarly  be  written  as 

=  7^-1  eM  (A I  -  A2Myl  d\  . 

2m  JT> 

Multiplying  by  V2Xf  and  obtaining  the  norm  of  both  sides,  we  have  the  estimate 

<  |^|  l  i'j,i  |  (v  -  r 1  i<«i  ■ 

Thus,  if  we  find  Mo  and  Ne  such  that  whenever  Re  A  >  v  +  e  and  N  >  Ne, 

(A  I-A2NYlVN  <M0, 

£(H) 

then  we  establish  (4.30).  This  uniform  boundedness  will  be  shown  by  decomposing  the 
resolvent  (A I  —  A2Xf)  into  its  components.  To  this  end,  consider  the  resolvent  equation 

where  (/, g)  €  R  and  {z^ ,z^)  €  Vixr .  Using  the  components  of  A  given  in  (4.19),  and 
letting  V2Af(f,g)  =  (fM ,g^),  we  obtain 

fAr  =  A  z?-z* 

g M  =  A jtf  +  AgzF  +  Afzf. 
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Eliminating  in  the  system  of  equations  above  and  using  the  assumption  A2  =  y*4i  +  A2 
we  have 


1  + 


A2 


dA/A-i 


Ay  + 

if*  ,  _ 

Ay  +  1  Ay  + 


+ 


Ay  + 

(O-1 


T«) 


ZM  - 
zl  — 


(4.32) 


r  (g*  +  Xf*  +  A?f*)  • 


To  simplify  notation,  let 


G  =  I  + 


Ay  + 


T^1 


1  +  j^TT  A-'AeWV) 


and 


F  =  7^~  +  1 

Ay  +  1  Ay  +  1 


(9  +  +  A  /) 


e  V  , 


and  denote  the  corresponding  finite  dimensional  expressions  by  G ^  and  respectively.  If 
we  show  that  ||y  is  bounded,  then  (4.30)  is  satisfied  where  is  given  by  G'K! z^  —  ■, 

for  ||(/,  <?)||-h  <  1,  Re  A  >  u  +  e  and  M  >  JVe.  The  next  step  is  then  to  show  that  the  inverse 
of  G  exists  and  is  bounded  in  £(V,V)  whenever  Re  A  >  v.  Since  (XI  —  A)  1  exists  and 
is  bounded  in  £(R,R)  whenever  Re  A  >  v,  then  for  every  (/,  g)  £  R,  we  can  solve  for 
(21,22)  €  R  such  that 


(A I  -  A) 


Zl 

Z2 


Solving  for  21  yields 


Zl 


1  + 


A 


Ay  +  1 


■A\  (A  H-  A2) 


if  ,  A 

Ay  +  1  Ay  + 


y  (9  +  A/  +  A2 /) 


and  ||2i||y  <  Mi  whenever  ||(/,  5)||^  <  1.  Thus,  G_1  exists  and  is  bounded  in  £(V.  V) 
whenever  Re  A  >  v.  Now,  we  consider  the  finite  dimensional  operators  and  show  bounded¬ 
ness  of  1 1  zJi  ||  y.  Note  first  that  G^ z^  can  be  expressed  as 


G*z?  = 


1  + 


AT\—1 


Ay  + 


T«) 


+ 


T«) 


Ay  + 


ZM 

zl 


=  1  + 


A2 


Ay  + 1 
A2 


Ai 1  + 


A 


7  A  A  zi 


Ay  + 
A 


Ay  +  1 

yMrMO-1)^ 


Ay  + 


Y  (a^A2  -  (A?)-1  A* )  zf. 


Since  G^ z^  =  F,M ,  we  have 

^  -  y^rr  (a/1  -  uf  r1)  -  YY^Y  (ay1^  -  (Afr'jg)  . 

Equivalently, 

^  =  G-1  (AI1  -  (A? r1) 

+  Ay  +  1  (-4l~1^2  ““  )_1Af )  |  • 
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Taking  the  V  norm  of  both  sides  and  using  Lemma  4.2  and  Lemma  4.3,  we  obtain 


w^\\v  <  I|G-1||£(v,v)(||^IIv  + 

{A^A2-P(f1AT1A2) ||£(W|K||v)  • 

From  the  convergences  in  Lemma  4.4  and  Lemma  4.5,  it  can  be  seen  that  pi^||y  is  bounded 
and  this  completes  the  proof.  ■ 

Remarks: 

1.  For  the  motivating  example  (2.4),  the  assumption  <72  =  707  +  a2  in  Lemma  4.1  is  sat¬ 
isfied  if  we  ignore  passive  patch  contributions  to  the  system  and  assume  that  the  den¬ 
sity,  Young’s  modulus,  air  damping  and  Kelvin- Voigt  damping  are  constant.  In  this  case, 
ai((p,ip)  =  fg  E I(f>" ip" dy,  a2(<p,  ip)  =  fg  Cnl(p"ip"dy  +  ca  fg  (pipdy  and  thus  7  =  Cd/E  and 
02  =  c0  fg  (pipdy,  where  02  satisfies 

Rea2(<p,<p)  =  ^UWh  >  ~\UWv  ■ 

The  operator  A2  generated  by  is  A2  =  caI  which  is  clearly  bounded.  We  then  write 
Ai 1  as  a  compact  operator  on  V  by  1  =  Af  li*i,  since  i  :  V  H  and  i*  :  H  V*  are 
compact  and  A^1  is  bounded  (due  to  (H2)).  Since  A2  is  bounded,  it  follows  that  A^1  A2 
is  compact  on  V. 

2.  If  passive  contributions  are  incorporated  in  the  model,  i.e.,  a±  and  a2  are  of  the  form  (2.10), 
Lemma  4.1  is  applicable  only  if  the  actuators  are  employed  in  such  a  manner  that  the  sum  of 
the  structure  and  actuator  stiffness  coefficients  is  a  multiple  of  the  sum  of  the  structure  and 
actuator  Kelvin- Voigt  damping  (this  is  necessary  to  satisfy  the  assumption  a2  =  701  +02). 
One  possible  way  to  achieve  this  is  to  embed  the  actuators  so  that  material  properties  of 
the  composite  structure  remain  constant  over  the  region  covered  by  the  patches. 

3.  The  assumption  of  exponential  stability  of  (A,B)  in  Theorem  4.2  is  guaranteed  by  Theo¬ 
rem  7.15  in  [2]. 

5  Concluding  Remarks 

We  presented  and  proved  Lemma  4.1  which  can  be  used  to  establish  uniform  stabilizability  of 
approximate  finite  dimensional  systems  arising  in  structural  systems.  In  the  context  presented, 
uniform  stabilizability  is  a  sufficient  condition  for  the  convergence  of  the  Riccati  solution  and 
control  gains.  As  demonstrated  in  the  remarks  above  and  in  [2,  Chapter  7],  the  conditions  in  the 
lemma  are  easily  verifiable  for  the  motivating  example  (2.4)  in  which  passive  patch  contributions 
are  ignored.  We  point  out  that  these  ideas  can  be  readily  extended  to  systems  involving  shells 
or  plates. 

For  models  incorporating  general  passive  patch  contributions,  the  assumption  a2  =  7 ay  +  a2 
with  <72  satisfying  the  hypothesis  in  Lemma  4.1  may  not  be  satisfied  and  the  authors  are  currently 
extending  the  lemma  to  include  this  more  general  case.  Such  an  extension  would  complete  the 
theory  for  convergence  of  Riccati  and  optimal  control  solutions  for  systems  with  no  exogenous 
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force  and  bounded  observation  operator  C.  Numerical  results  demonstrating  convergence  of  sub- 
optimal  controls  for  thin  shell  structures  incorporating  passive  piezoceramic  patch  contributions 
have  been  demonstrated  in  [12,  4,  5]. 
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